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Abstract 

In this paper a new block-structure method is presented for the solution of 
the well-known from gravity theory matrix system of equations g%jg^ k = 5 k (with 
respect to the unknown covariant components gij and by known contravariant ones 
gi ) by transforming this matrix system into a linear algebraic system of equations 
in the general N— dimensional case. Although powerful computer methods exist 
for the solution of this problem for a given (fixed) dimension of the matrices g lJ 
and especially for numerical elements of g %3 , the structure of the obtained linear 
algebraic system in the general N— dimensional case and for arbitrary elements of 
g %3 (functions) has not been known. 

The proposed new analytical block-structure method for the case of symmetrical 
matrices gij and g jk (the standard case in gravity theory) is based on the construc- 
tion of a block-structure matrix, whose "elements" are again matrices. The method 
allows to obtain the structure of this linear system in the general N— dimensional 
case, after multiplication (to the left) with the transponed matrix. 

Some arguments are given why the proposed method may be applied, after some 
refinement and generalization for the case of non-symmetrical matrices gij and g jk , 
for finding the graviton modes in the Kaluza-Klein expansion in theories with extra 
dimensions. 



*Electronic mail: bogdan@theor.jinr.ru 
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1 INTRODUCTION 



The system of equations g^g-** = 5 k is well-known and important in gravity theory, since 
from it the contravariant metric tensor components are determined under known covariant 
ones. In spite of its importance, this system is not well-understood as a mathematical 
object. For example, in Bergman's book on General Relativity [1] (Ch. 5, eq. 5.64) 
it is written that " if the determinant gij is not equal to zero, then a multitude of new 
variables g^ k can be introduced according to the relation gifg^ k = 5 k ". No doubt, similar 
statements can be seen in other books. Bergmann's statement would have been true if 
this system is a trivial linear algebraic system of equations with a number of equations 
equal to the number of variables. But in fact, this is a predetermined system of equations 



with a number of equations n 2 greater than the number ^ ^ j + n of unknown variables. 

This is the first reason why Bergmann' s statement does not hold. Secondly, this is not a 
linear algebraic system of equations, but an operator (matrix) system of equations [2] of 
the type AX = B, where X is no longer a vector-column, but a matrix. Consequently, 
the linear algebra theorems do not hold any more. 

One of the basic new results in the present paper is that the initially given operator 
system of equations can be reduced to a linear algebraic system of equations BY = T, 
where the unknown elements Y represent the matrix elements X, arranged in the form of 
a vector-column. At first, the system will be written in the form AY = 5 k , where A is a 

2 j + n xn 2 matrix. Further, non-quadratic matrices will be denoted 

as a x b, where a will be the number of rows in the matrix and b will be the number of 
columns. 

In the conclusion, a possible application of the proposed mathematical method will 
be briefly commented in reference to the problem about finding the graviton's modes 
in the Kaluza-Klein's expansion under toroidal compactification in theories with extra 
dimensions. 



non-quadratic 



2 Some general properties of the system gi^ k = h\ 
in the general n— dimensional case. 

Since the number of equations is greater than the number of variables, from the n 2 equa- 
tions one can select n 2 — n equations with different values of i and k, for which the 
right-hand side (R. H. S.) will be zero. Yet the number of the chosen equations remains 

to be greater than the number of variables ( ^ ) + n, which is confirmed by the equality 



(n 2 — n) — 



+ n 



n{n — 3) 



(2.1) 
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fulfilled for n > 3. Only for the case n — 3 , the number of the equations n 2 — n = 6 with 

n 



a zero R. H. S. becomes exactly equal to the number of variables 
^- = 6. Therefore for the case of an arbitrary n from these n 2 



choose 



+ n 



equations. In each row of the 



+ n 



x 



+ n 



matrix 



of coefficient functions of this system there will be only n functions g jk , the rest 



n 



of the elements will be zero. The determinant of the matrix will be equal to zero, a proof 
of which in the general case of an arbitrary n (and also for n = 3) will be presented in the 

system of equations with a zero determinant are arbitrary. 
Now we are left with w( "~ 3 - ) equations (again with ^ ^ 

R. H. S. (we shall call it the first system of equations) and with another n equations (the 
second system) with a R. H. S. of each equation, equal to 1. Since the number of variables 



n dimensional homogeneous 
+ n variables) with a zero 



n 



+ n in the first system is greater than the number n ^ r \ — of equations and 



+ n- 



n(n 



3) 



2n 



(2.2) 



again one can treat as unknown only - variables and transfer the rest 2n variables in 



the R. H. S., which will become different from zero. Analogously, for the second system 

n> \ n(n— 1) 

2 



one can treat as unknown only n variables and transfer the rest 



variables 



in the R. H. S. If this R. H. S. is different from zero and moreover, the determinant of the 
coefficient functions <f J is also different from zero, then one can find unique solutions for 
these n variables g^. Now comes the most important point of the proof: Since 



n(n — 1) 
n + — > 2n 



(2.3) 



one can take all n of the uniquely found solutions plus n more unfixed (freely varied) 
variables from the R. H.S. of the second system and "plunge" them into the R. H. S. of 
the first system. Let us remember that the first system has " ( -"~ 3 ' > unknown variables, 

but one may note that ^™^j+n — 2n= n ( n ~ 3 ^ variables from the second system have 

not been transferred in the R. H. S. of the first system. Therefore, one can choose these 
n( - n ~ 3 - ) variables to be the unknown variables for the first system, and if the determinant 
of the coefficient functions is nonzero, an unique solution can be found for them. 

As a whole , one would have a maximum of n+ w(w ~ 3 - > = " ( -"~ 1 - > uniquely fixed variables 
and the rest of the variables may be varied freely. 
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3 BLOCK STRUCTURE METHOD FOR THE 
PARTIAL n = 3 CASE 



The purpose of the present section will be to develop a method for solving the system of 
equations g^g^ = 5* for the case n — 3 . Since in principle the system for n = 3 can 
be solved in an elementary manner, the aim will be not to find another more convenient 
method, but rather than that find a method, which can further be generalized to higher 
dimensions. 

In matrix notations for the unknown variables gn, #22, #33 , <?i2, #23 and #13 the system 
gijg jk — $1 of nine equations (for the different indices i and k) can be written as 

AX = E , (3.1) 

where X T further will denote the transponed vector 

X T = (gn, #22, #33,312, #23, #13) • (3.2) 

and E T - the transponed 9— vector: 

E T = (1,0, 0,0, 1,1, 0,0,1) . (3.3) 

The 6x9 matrix A has the following interesting block structure: 



(3.4) 



where (s = 1, 2, 3) the matrices P s and Q s are the following: 

gsl g s2 g s3 \ / 

P s = I g sl ; Qs=\ g 2s g 3s \ (3.5) 

g sl J \ g 2s g 3s 

In order to find the solution X = A~ l E oi the system (3.1), one has to find the inverse 
matrix A^ 1 . For the case of quadratic matrices with the block structure 

m^(cd) • < 36 > 

where A,B,C and D are n x n, q x n,n x q and q x q matrices correspondingly, the so 
called Frobenius formulae [2] for finding the inverse matrix M~ l is valid 

M -x _[A- 1 + A-'BH-'CA- 1 -A-'BH- 1 \ 
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where H is the matrix: 



H = D -CA~ X B . 



(3.8) 



In the present case, the Frobenius formulae cannot be applied to the block matrix (3.4), 
since it is not a quadratic one. However, if X is a solution of the system (3.1), then it is 
a solution also of the equation A 1 AX = EA T , where the 6x9 matrix A, multiplied to 
the left with its transponed one, gives already the quadratic 9x9 matrix A T A. Further 
it shall be demonstrated how the Frobenius inversion formulae can be applied twice in 
respect to A 7 ""A 

The matrix A T A can be calculated to be the following block matrix: 



A A 



\ 



i=i 

3 



3=1 
3 



* E Ql p k EQlQi i 
\ k=i i=i / 



(3.9) 



and Pi,Qj are the corresponding matrices (3.5) and their transponed ones. The block 
matrices in (3.9) are found to be the following 3x3 matrices, which shall further be 
identified with the corresponding block - matrices A, B, C and D in (3.6): 



i=l 



Jl„i2 



/ 3 3 

' E^ 1 ) 2 E 

i=l i=l 

E^V 2 E[(^) 2 + (^ 

i=l i=l 
3 3 

\ i=l i=l 



E g a g* 3 
i=i 

E 9 l2 9 13 
i=i 

Jl\2 



\ 



E t^ 1 ) 2 + (9 t3 ) 

i=l 



(3.10) 



3=1 3=1 





g2j g 3j 

tf 2 g^ 



k=l 



k=l 




1=1 



2i r ,n 



i=i 



o 



( Z(g 21 ) 2 
i=i 

Eg 2l g 31 E [(g 21 ) 2 + (g 31 ) 2 } Eg 2l g 

i=i i=i i=i 

o EsV E(// :i/ ) 2 

1=1 i=l 



3/ 



/ 



(3.11) 



(3.12) 



(3.13) 



Note that the diagonal block - matrices A and D have non - zero determinants (even 
if g l i = dX l dX^) ) while the non - diagonal block - matrices B and D have zero - 
determinants. However, in order to apply the Frobenius formulae for inverting the matrix 
(3.9) it is sufficient to have as invertible only the matrix A (and of course the matrix H). 
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4 MODIFICATION OF THE BLOCK STRUCTURE 
OF THE MATRIX A- THE n = 3 CASE AND THE 
GENERAL TV— DIMENSIONAL CASE 

The above presented method has nevertheless the following shortcomings: 

1. It deals with an rectangular p x q matrix A for a system of equations with p = 

ft \ 

+ n variables and q = n 2 equations. At the same time it would have been much 



2 

better to deal with a quadratic matrix at the beginning. 

2. The block - matrix A contains two types of matrices P s and Q s , while it would be 
more convenient to have just one type of an elementary "constituent" with a definite 
structure, where the indice i denotes the corresponding column in the block matrix A 
(i.e. the number of the column, containing block - matrices) and the indice k - the 
corresponding row of block - matrices. 

3. The matrix Q s , given by the formulae (3.4) has a zero determinant. 

To avoid these shortcomings, let us define an extended 9— dimensional vector Y, whose 
transponed one is 

Y T = (911,913, 913,921,922, 923, 931,932, 933) , (4.1) 

where the elements 921,932 and g$i formally shall be considered unknown, although 
they are equal to their symmetric counterparts gi2,923 and gi 3 . Then the system of equa- 
tions can be written as 

MY = 1 , (4.2) 

— T 

where 1 is the transponed 9 -dimensional vector 

1 T = (1,0, 0,0, 1,0, 0,0,1) (4.3) 
and the 9x9 matrix M has the following block structure 



M 



( E[ 2) E[ 3) \ 

Ef ] E® 
\ E U E P Ef J 



(4.4) 



The elementary 3x3 block matrices E^\ E^ and Ej® (k = 1,2,3 denotes the number 
of the row) in each column are the following 

E k^\\ g kl | , (4.5) 



2g kl 


g k2 


g k3 





g kl 











gkl 


g 2k 








g lk 


2g 2k 


g 3k 








g 2k 



Ef ] = \\ <j u 2,;* </* | (l.G) 



6 



k 2 



9 



3A- 







5 







,2fc 



(4.7) 

gM g ^ 

From (4.5 - 4.7) it is seen that depending on the indice k of the elementary matrices row 



and of the indice (s) for the elementary matrices column, the matrix E k 
dimensional case can be written as 



for the n— 



E 



(n-1) 



/ # fcs 

# fcs 0.. 

V 








o \ 





(n-l)g 



ks 



-,kn 



(4.8) 



/ 



In other words, on the main diagonal of the matrix E^ the elements g ks are situated, 

on the s— th row - the elements ( g kl , g k2 , . .. (n — l)g ks , g kn ) with an element 

(n — l)g ks on the s— th row and s— th column. This block structure (with some slight modi- 
fications) shall be obtained also for the n— dimensional case, and thus the 3— dimensional 
case really helps to make the corresponding generalization for the n— dimensional case. 

Another advantage of the block - matrix representation is that it gives a possibility 
to apply twice the Frobenius formulae. Correspondingly, in the n— dimensional case the 
Frobenius formulae will be applied (n — 1) times. 



5 BLOCK MATRIX STRUCTURE IN THE N- 
DIMENSIONAL CASE 



Following the same algorithm as in the preceeding subsections, we shall try to find the 
block structure of the system of equations gijg jk = 5 k in the n— dimensional case. For 

/ \ 



% — k the system can be written as AXk 



vector X k is 



and the matrix A is 



1 

V - ) 



, where " 1" is on the A;— th place, the 



Xk = (<7fci> <7fc2, ,9kn) 



( g 11 g 12 g ln \ 



A = 



g 21 g 22 



nl „n2 



Jin 



\ 9 nl 9 



9 



(5.1) 



(5.2) 
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The corresponding vectors Xi,X 2 , ,X k represent the rows of the symmetric matrix 

N 

( 011 #12 gin \ 

921 922 92n 



N = 



(5.3) 



\ 9nl 9n2 •••• 9nn ) 

in which the unknown variables are in the lower triangular (half) part of the matrix 
(denoted by N tr ). 



Let us construct a 



x n 2 dimensional matrix B, which will multiply a w(w 2 +1 - > 
dimensional vector Y, formed by joining all the consequent rows of the triangular matrix 

N tr 

Y = (gn,gi2, , gin, 922, #23, , g2n, ,g n i,g n 2, ■■■ ■•; ynn ) • (5-4) 

Correspondingly the matrix B will have the following block triangular structure: 



B = 



( B n 

B 2 \ B 2 2 



ni 



B 



n2 



\ 



B nn J 



(5.5) 



where each of the block matrices B kk on the diagonal is an n x (n — k + 1) (i.e. (n — k + 1) 
columns and n rows) dimensional matrix, obtained from the matrix A by removing the 
first (k — 1) columns. For example, B n = A, but 



B> 



22 



/ g 12 g 13 

g 22 g 23 



y g n2 g n3 



The corresponding matrix B kk will be 



B kk = 



/ g lk g l,(k+l) 
g 2k 2,(fc+l) 



g m \ 

g 2n 
g nn J 



g l,n \ 
g 2,n 



(5.6) 



\ g n ^ ^ ri '(^+^) 



(5.7) 



g nn j 



The block - matrices B sk (s > k) are n x (n — fc + 1) dimensional ones with just one 
nonzero column (the s — k + 1 column) with the elements g kl , g k2 , 



B s k = 



/ o 




V o 



/ti 

,fc2 



9 



9 
9" 

fc,(n-l) 



9 



... \ 

.... 

.... 

.... J 



(5.8) 
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Let us now have a look at the complex of neighbouring block matrices around the main 
block diagonal 

K (k) = ( Bfc-i,*-i B k-i,k \ 

\ £>k,k-l B kk J 

For illustration of the block matrix multiplication and in order to derive some useful 
formulaes, let us calculate K^ T K^ k \ which will be equal to 



B p-l,p-l B k-l,k-l + B p ,p-l B k,k-l B p,p-l B k,k 



B pp B k,k-l 



BLB 



pp 



kk 



The corresponding terms in the above matrix are: 

B p-l,p-l B k-l,k-l 



g 

g l,P 



\ 9 hn 



2,(p-l) 


gU,{p-l) 


9 2 > p 




g 2,n 


gUn 


Pu 


Pl2 



g l,(k-l) l,fc 

g 2,(k-l) g 2,k 



,n,fc— 1 



V 9"'" ' g 

Pn-k+2 



n,k 



g l,n \ 
g 2,n 

cj nn j 



(5.10) 



(5.11) 



\ Pn-p+2,1 Pn-p+2,2 Pn-p+2,n-k+2 ) 

where all the elements of the matrix are nonzero and the element P sr on the s— th row 

n 

and on the r— th column is equal to P sr = Yl g % ^ +s ~ 2 g l ^ k+r ^ 2 , 

i=i 



/ 



B p,p-l B k,k-l 







giP-V'ig^- 1 )'' 1 



\ 




1=1 



\ 







(5.12) 



/ 



where x and B^k-i are (n — p + 2) x n and n x (n — A; + 2) matrices and the 

n 

resulting n x n matrix has only one nonzero element Yl g^ -1 ^^ -1 )' 1 on the second row 

i=\ 

and on the second column. Next let us calculate the matrix which is a product 

of the (n—p+1) x n matrix B pp and the n x (n — + 2) matrix B kk _i with the only 
nonzero second column: 



B pp B Kk _i 



( o 





F\2 

F22 



... 0\ 
... 



(5.13) 



V F„_ p+li2 ... J 
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The resulting matrix B pp B kk _i has a dimension (n — p + 1) x {n — k + 2) with the 



only nonzero second column with an element on the r— th row and on the 2— nd column 

n 

F r2 = £ gi,(p+r-l) g k-l,i 



i=l 



Next let us find the matrix which is a product of the (n — p + 2) x n 

matrix B^ v _ x and the n x (n — k + 1) matrix B k> k : 



B P,P-i B k,k 



( 



E 



V 










n 






\ 

n 

g < J>- 1 )^ i g i ' n 





o ) 



(5.14) 



where in the last formulae we have used that the obtained matrix has (n — k + 1) columns 

n 

and therefore the indice k + s — 1 in the expression for the element ^ gip-tf^ghik+s-i) on 

i=i 

the 2— nd row and s— th column ranges from k to k + s — 1 — k + (n — k + 1) — 1 = n. 

It remains only to calculate the matrix B pp Bkk, but it is the same as (5.11), this time 
with an element 



,p+s-lgi,k+r— 1 



(5.15) 



?=i 



on the s— th row and on the r— th column. 

Using the above developed techniques for matrix multiplication, let us calculate the 
^ x ^±i) matrix £? T £? (recall - B is an n 2 x ^±i! matrix and S T is an ^±i! x n 2 ma- 
trix), which is the n— dimensional analogue of the matrix (3.8). Taking into consideration 
(5.9), one has 



/ n n n \ 

( J2 B Tl B H J2 B Tl B i2 ••• Y. B Il B in 



B B 



i=l 



i=2 
n 



t=n 
n 



E B i2 B H E B i2 B 



i2 



i=2 



\ E B Tn B H 
\ i=n 



The above matrix contains three types of elements: 



E B i2 B n 



E B In B v. 



(5.16) 



1-st type. Elements below the block diagonal of the type B^-B^ with k < j 



a=] 



Carrying out the matrix multiplication and for the moment not taking the summation 



10 



over a, we find 



( 

.. 

.. 
\ 

..ji „ki 



E 9 ji 9 ki 

i=l 



\ 







.0 

J 



(5.17) 



where the only nonzero element E gi l g kl in the matrix is on the (a — j + 1) row and on 



the (a — /c + 1) column and the matrices and £? a £ are of dimensions (n — j + 1) x n 



and n x (n — k + 1) correspondingly. 

n 

Since the first term in the sum E B^.B ak will be B^ a B ak , let us find it, performing 

a=jf 

the same kind of matrix multiplication as in (5.13): 



bLb, 



ak 







E 

i=l 



Act , t ki 



o .. Es 

i=i 

o .. 

n 



,a+r—lgk,; 



i,n„k,i 



o o Ert 
\ i=i 





.. 

.. .0 

.. 



\ 



(5.18) 



where the only nonzero column is the (a — k + 1) one and the element in this column 

n 

and on the r— th row is E g' t ' a+r ^ 1 g k '\ The matrices B^ a and B ak are of dimensions 



i=i 



(n — a + 1) x n and n x (n — /c + 1) correspondingly and the resulting matrix B^ a B ak is 
(n — a + 1) x (n — A; + 1) dimensional. 

2-nd type. Elements above the block diagonal of the type E B^.B ak with k > j [The 

a=k 

summation indice a takes at first the value of that indice (k or j), which is greater than 
the other]. The first term in the above sum is B^B aa , which can be found to be 



B aj B aa 



\ 



^2 gj> i g 1 '( a + i - 1 ) ... ^2 gjjgPXa+i- 1 ) .. ^2 g3,ign,(a+i-l) 
i=l i=l i=l 



\ 



















(5.19) 



/ 



The above (n — j + 1) x {n — a + 1) matrix B^.B aa contains a nonzero (a — j + 1) row 
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with an element in the p— th column, equal to 

n 

i=i 

3-rd type. Elements situated on the block diagonal of the type 

n n 
^2 B ak B ak = Bl a B aa + ^ B lk B o 



(5.20) 



>ak 



(5.21) 



a=k 



a=k+l 



Similarly to the calculation of (5.15), the first term in (5.21) (B^ a is an (n — a + 1) x n 
matrix ; B aa is an nx (n—a+1) matrix) is found to be the following x (n— a+1) 

matrix (also, a — k) : 



B Ia B »a 



( N u 

N 21 



N 12 
N 22 



Nl,(n-a+l) \ 



2,(n-a+l) 



+1,1 



N, 



(5.22) 



n— a+l,n— a+1 



where 



AT, 



P9 



,a+p-l^i,a+<j-l 



(5.23) 



i=l 



The second term in (5.21) is in fact the (n — k + 1) x (n — k + 1) matrix (5.17) for j = k. 
The summation over a from a = k + 1 to n will give a diagonal (n — k + 1) x (n — fc + 1) 

n 

matrix with an element G a -fc+i, Q -fc+i = ^2(9 k ' % ) 2 011 the (a — k + 1) row and on the 

(a — k + 1) column, which will be situated on the main (block) diagonal from a — k + 1 
to n. Since for a = k + 1 we have a — k J r l = k + l — k + 1 = 2 and for a = n we 

n 

have also a — fc + 1 — n — k + 1, this means that the matrix E B^ k B ak will have the 

a=k+l 

following structure: 



/ 

n 

E(/'T 



X] B ak B ak = 



a=k+l 
















E(^) 2 










\ 



\ 



(5.24) 



.... Z(9 k 'T / 

i=i / 

Therefore, summing up the matrices (5.22) and (5.25), one obtains the general structure 
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of the matrices (5.21) Yl B^ k B ak on the block diagonal: 

a=k 



( 



^2 B ak B ak 



a=k 



N 2 i 



N 12 
N22 + G22 



^2,(n-a+l) 



\ 



+1,1 



— a+l,n— a+1 + G22 J 



(5.25) 



where N pq and G22 are given by expressions (5.23) and G*2,2 = E(3 ) 2 respectively. 



Consequently (for p > 2) 

A^p + G22 



i=i 



Eb 



i,(fc+p-l) i,(fc+p-l) 



(5.26) 



Let us find now the structure of the off - diagonal block matrices, situated below the 
diagonal in the block matrix (5.16). Each such an element can be decomposed as 



^2 B aj B ak = B aa B ak + ^ B ^3 Bak ' 
a=j(j>k) a=j+l 



(5.27) 



The first term in (5.27) is the already known (n — a + 1) x {n — a + 1) matrix (5.18) with 
the only nonzero (j — k + 1) column. The second term is the sum from a = j + l to a = n 
of the (n — a + 1) x (n — k+1) matrices B^jB ak , already calculated in (5.17) and having 

n 

the only nonzero element gi l g kl on the (a — j + 1) row and on the (a — k + 1) column. 

i=i 

n 

The summation over a from a = j + ltoa = n means that in the sum YJ B^B ak the 

a=j+l 

n 

element YJ g^ l g kt will appear beginning from the row a—j + 1 = 2 (for a = j + 1) up to the 

i=i 

row a — j + 1 = n — j + 1 (for a — n), which in fact is the last row. Correspondingly, the 
same element will appear beginning from the column a — k + 1 = j — k + 2 (for a — j + 1) 
and ending at the column a — k + 1 = n — k + 1 (at a = n), which is the last column. In 
other words, the summation of the matrices B^B ak , containing one element, effectively 
results in a matrix, filled up from the second row to the end and from the (j — k + 2) 
column to the end: 



n 



n 



/ , BljB ak = ^ 



a=j+l 



a=j+l 



( 



\ 







<£gjigki 
1=1 





\ 



J 
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/ 





column (j — k + 2)... 

n 

i=i 



o 

n 

E g ji g ki 

i=l 



\ 



Eg Jl g kl 
i=i 

n 

E ^V* 
i=i 



^gjigki 
i=l 

•• E^V* 
i=i 



(5.28) 



/ 



Now recall that the matrix B^ a B ak (5.18) had a nonzero (j — k + 1) column, so therefore 

n 

the structure of the whole matrix E B?.B ak (j > k) below the diagonal is similar to 

a=j 

(5.28), but with the additional (j — k + 1) column: 

( 





0. 



V 



0. 



n 

i=l 
n 

^2 gi,a+lgk,i 
i=l 








n 

E^V 1 • 
i=i 


• T,g ji 
i=i 


n 

gi,a+r—l gk,i 

i=l 


E • 

i=i 


n 

• E^ 

i=i 


n 

i=i 


E^V' • 
i=i 


n 

• E^ 

i=i 



(5.29) 



/ 



In a completely analogous way the elements above the block diagonal in (5.16) can 

n 

be found. These elements E Bj-B ak (k > j) can be decomposed as 



a=k 



^2 B IjB a k = B kj B kk + ^2 B T aj B ak 

a=k(k>j) a=k+l 



(5.30) 



This formulae is similar to (5.27) for the below - diagonal elements, but here in (5.30) we 
have the matrix B^B kk instead of the matrix B^ a B ak and the summation over the indice 
a in the second term is from a = k + 1 to a = n instead of a = j + 1 to a = n in (5.27). 
The first term in (5.30) B^B kk is the (n-j + 1) x (ra-fc + 1) matrix (5.19) (with a = k) 
with the only nonzero (k — j + 1) row with the elements 



XyV' ^ 5 ••••>>V- / ,; ' : E^V'^ (5.31) 

i=i i=i i=i 

The second term in (5.30) is again the sum of the matrices (5.17). The only nonzero 

n 

element E g^g kt will now appear in the final sum from the a — j + 1 = k — j + 2 row 

i=i 

(a — k + 1) until the k — j + 1 = n — j + 1 (a — n) row, which is the last one. Also, 
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the same element will fill up the columns from a — k + 1 = 2 (a — k + 1) to the column 
a — k+1 = n — k+1 (a = n), which is also the last one. Therefore, summing up the two 
terms in (5.30), one obtains the following matrix for the above - diagonal block terms, in 

n 

which the (k — j + 1) row is filled up with the elements E gi^gf'^ 1 ^ and from the next 

i=i 

(k — j + 2) row to the end and from the second column to the end column the matrix is 

n 

filled up with the other element E g il g k% : 

i=i 



/ 



E gi' i g 1 '( a+i ~ 1 ') 
i=i 





0. 



1=1 

n 

E 9 ji 9 h 

i=l 



^ki 



E 9 ji 9 ki 

i=l 



0.. 



E g^ i g p ^ a+i ~ l ^> 
i=i 



n 

E 9 i,n 9 

i=l 



E ^ w 

i=l 



E ^ w 
i=i 



gj,ign,(a+i-l) 

i=l 

E # j V* 
i=i 



E fc * 
i=i 

(5.32) 



Let us now summarize the obtained results for the n— dimensional case. The (prede- 
termined) system of equations g^ k = S k was represented as BY = 1 , where 1 is an n 2 
dimensional vector, whose transponed one is defined as 1 



'l T \l T2 



, l Tn ) and 

l Tn are defined 



the corresponding transponed n— dimensional vectors 1 ,1 , 

2^1 T2 

as follows: 1 = (1,0, ,0),1 = (0,1,0, 0) and the k— th transponed vector 



-rTk 



(0,0, , 1,0, ...0) contains the number 1 on the k— th place. The n 2 x 



n(n+l) 



ma- 



trix B in terms of the elementary block matrices By and the w( -" 2 +1 ' ) dimensional vector Y 
were defined by formulaes (5.5 - 5.8) and (5.4) respectively. In order to solve the system, 
we multiplied it to the left with the transponed matrix B T and thus the solution for the 
vector Y in matrix notations can be found as Y = (B T B)~ 1 B T 1, where the expression 
for B T 1 can easily be found to be 





( £n 


#21 ••• 


... Bl x 






\ 


B T 1 = 





B22 ■■ 


■■ B% 2 




I 2 











■ B T nn 


) 




/ 



( 



E BIT 

i=i 

n 

E BIT 

i=i 

n 

E BJJ 

1=1 



\ 



/ 



(5.33) 



The expressions for the elements of the vector in the R. H. S. can also be found, but this 
will not be performed here and will be left for the interested reader. Let us remind again 
that each "element" of this vector is in itself an n— dimensional vector and m = w(w 2 +1 ' > 
gives the number of the "block" elements of this vector - column. The number m should 
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be an integer number, but this requirement can be fulfilled and this will be commented 
in the next subsection. 

The main result of this subsection are contained in the expressions (5.25), (5.29) and 
(5.32) for the elements of the matrix B T B on the block diagonal, below the block diagonal 
and above the block diagonal respectively. In such a way the detailed structure of the 
matrtix B T B is known in terms of the elementary constituent block matrices B^. The 
structure of the matrix B T B is important for the following reasons: 

1. If one chooses the contravariant metric tensor in the form of a factorized product 
<f J = dX x dX\ then one can answer the question whether it is possible the rows of the 
matrix B T B to have some common multiplier. But from the above expressions and since 

n 

in each "cell" of the elementary block matrices one has a summation of the kind g l,n g k ' 1 , 

i=i 

it is evident that such a common multiplier does not exist and therefore the rank of the 
matrix B T B cannot be lowered due to the above made choice of g' lj . 

2. The quadratic structure of the matrix B T B gives an opportunity to apply the 
Frobenius formulae for finding the inverse matrix (B T B)~ l ,as previously discussed. 

6 BLOCK MATRIX REPRESENTATION OF THE 

HOMOGENEOUS SYSTEM OF EQUATIONS 
WITH A ZERO R. H. S. 

Earlier it was shown that for the case n = 3 the sub - system of equations g^g^ = with 
a zero R. H. S. has a determinant of coefficient (functions), equal to zero. The question 
which naturally arises is whether this property is valid only for the n = 3 case and is it 
valid for the n— dimensional case? 

Below it shall be proved that this can be done for the general case. Namely, it shall 
be established that from the n 2 — n equations with a zero R. H. S. there may be chosen 

a sub-system of ^ ^ ^ + n = " ( " 2 +1 ' > equations with a determinant of coefficients, equal 

to zero. It shall be stressed that the proof will be that such a system exists (i.e. can 
be chosen) and not that all other choices of the subsystem of equations will also satisfy 
this requirement. In fact, investigating under what other choices of the sub - system of 
equations with a zero R. H. S. this property will be preserved, represents an interesting 
problem for further research. 

For the purpose, let us again use the block matrix representation (5.5). Then, in order 
to obtain the system of equations with a zero R. H. S. , one has to remove the first row 
(for % — k — 1) from the system g nk gu = 5™ of the first n equations (i.e. for n = i), then 
the second row from the second system of n equations and so on, one has to remove the 
k— th row from the k— th system of n equations. Further, in order to receive again a block 
matrix structure with (elementary) submatrices with n rows, one should add the first row 
of the second system of n equations as the last row of the first system of equations. In 
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effect, since the first two rows of the second system have been removed, one should add 
the first two rows of the third system of n equations as the last two rows of the second 
system. Therefore, since also the third row in the third system ( for % — k — 3) has been 
removed, it has a total of three rows removed, and subsequently three last rows have to be 
added from the fourth system. Continuing in the same manner, from the k— th matrix Bkk 
on the block diagonal we have the first k— rows from B kk removed and also k - last rows 
added, which should be taken from the below - diagonal matrix B k+ i k . Since according 
to (5.8) the n x (n — k + 1) matrix B sk has a nonzero (s — k + 1) column, the matrix B k+ i :k 
will have a nonzero k + 1 — k + 1 — 2 column. Therefore the new transformed in this way 
n x (n — k + 1) matrix, which will be denoted as Bkk , will have the following structure: 



f 9 



B, 



kk 



k+l,k gk+l,k+l 
gk+2,k gk+2,k+l 



-,nk 



^n,(fc+l) 
„kl 



9 



kk 



gk+l,n 
gk+2,n 







(6.1) 



/ 



In the same way, the below - diagonal transformed matrix B k+1>k , obtained from B k+lk 
by removing its first k and adding k rows from B k+2<k , will be of the following kind: 



B 



k+l,k 








,k,n 







9 







\ 






(1) 



\ g K ^ ..00/ 

Since the block matrix (5.5) has a triangular structure, for our further purposes only the 
structure of the block - diagonal matrices Bkk will be relevent. 

Next, our goal will be to divide the the block - matrix (5.5) into elemetary block 
matrices with an equal number (n) of rows and columns. Let us remind once again that 
the block - matrix (5.5) contained elementary block - matrices with an unequal number of 

columns - n, (n — 1), (n — 2) correspondingly. For the purpose, we shall take one (left) 

column from the block matrix with (n — 2) columns and transfer it to the left to the block 
matrice with (n — 1) columns. As a result, the block matrices on the second block matrix 
column (B. M. C.) will already contain n columns. Since in the block matrix column one 
column has been transfered, one has to add 3 columns from the (n — 3)— rd block matrix 
column to the (n — 2)— nd block matrix column in order to obtain again a block matrix 
column, consisting of elementary n x n matrices. Continuing in the same manner with 
the (n — 3)— rd B. M. C, one has to add to its right end 6 columns from the (n — 4)— th 
B. M. C. Now let us write down the numbers of the corresponding block columns and 
below with a (— ) sign the number of columns, transferred to the (neighbouring) B. M. C. 
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; with a (+) sign the number of columns, joined to the B. M.C. (to the right side) from 
the neighbouring (right) B. M. C. 

(n-1) ; (n-2) ; (n-3) ; (n-4) ; (n-5) 

+1; -1 +3 ; -3 +6 ; -6 +10 ; -10 +15 

Now let us look at the numbers with a minus sign, which form the following number 
sequence (with the corresponding number in the sequence denoted): 



1,3,6,10,15,21 

1,2,3,4,5,6, (6.2) 

It is trivial to note that each number in the sequence (upper row) is in fact a sum 
of the corresponding numbers in the lower row up to that number. For example, the 
number 21 in the sequence (upper row) can be represented as a sum of the numbers in 
the sequence (lower row): 21 = 1 + 2 + 3 + 4 + 5 + 6. The same with the number 10 
==>- 10 = 1 + 2+3 + 4. Therefore, to the sequence number k in the low row will correspond 

the number fc ( fc+1 ) in the upper row, which is the sum of the first k numbers in the low 
row. Since the number k in the lower (6.2) corresponds to the (n — k) — th block column, 
the corresponding number will be k ^ k ~^ ; and it will correspond to the number of columns, 
which have to be removed from the diagonal block matrix B k+ i^ k+i y At the same time, 
to the right end one should add k + fc ^ fc ~ 1 ) = M*±I2 columns. This number is exactly equal 
to the number of left columns, removed from the (right) neighbouring matrix B^+i),(k+2)- 
This can serve also as a consistency check that the performed calculation is consistent 
and correct. 

From the matrix (6.1) for B k+1 ^ k+1 ^ we have to delete the first s + 1 left columns, the 

first (upper) elements of which begin with the elements g( k+2 )'( k+1 )- g(fc+2),(fc+2) ^ g(k+2), 

where s = k ^ k ~^ , Now let us denote by B^ + x),{k+i) the matrix B k+1 ^ k+1 ^ with M*±ll 
columns removed and k i k ^ll right columns added. The upper elements of the (left) re- 
maining columns will be g( k + 2 )^( s + k + 1 ) ^ g(k+2),( s +k+2) ^ ^(fc+2),^ wriere 



s + k + 1 



k(k- 1) 



+ k + l 



k(k+l) 



+ 1 



(6.3) 



2 2 

and the matrix B^+i),(k+i) wm contain n — k ( k +^ left nonzero columns. Therefore the 
n x n matrix will have the following structure: 

-B(fe+i),(fc+i) = ( lI ' ( 6 - 4 ) 



where the (n — k + 1) x 



n 



fc(fc+i) 

2 



matrix Li is 



/ £? (fc+2),(M*+i)+l ! 



(fc+3),(M^tl)+l) 



g (k+2),n \ 
(fe+3),n 



9 



k(k+i) 



+i 



9" 



(6.5) 
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and L 2 ,L 3 and L 4 are zero matrices of dimensions (n — k + 1) x 



fc(fc-i) 

2 



(fc + 1) X 



77. 



fc(fc+l) 
2 



and (/c 
Li and L3 contain 



+ 1) x 



fc(fc+i) 
2 



n 



fc(fc+i) 
2 



correspondingly. Note also that the block matrices 

columns of 



fc(fc+i) 
2 



nonzero columns, and the remaining 
the matrices L 2 and L4 are exactly equal to the number of zero columns, added to the 
right side of the matrix B^+i),(k+i) from the neighbouring matrix B^+i),(k+2) ■ Since this 
result depends on the initial structure of the matrix and on the expression 

(6.3) (which are both independent on the number of removed right columns), this also 
confirms the consistency of the calculation. 

Note that the block structure of the matrix (6.4) has been revealed on the base of the 
assumption that the elements ...gr( fe + 1 )>( fe + 1 ) i n the last (k+ 1) rows and the second 

column in the matrix (6.1) will be among the first removed to the left (and outside the 
matrix) columns. However, for s = (i.e. k — 1) the elements in the second column of 
the last two rows of the matrix B 2 ,2 will contain the elements g 21 ,g 22 - Therefore, after 
removing the first column in the matrix (6.1) Bfk+i),(k+i) (f°r k — 1) and adding to the 
right one zero column, the obtained structure of the nxn matrix -62,2 will be the following: 



-82,2 = 



( 9 33 


9 U ■ 


• 9 3n 





\ 




9 U ■ 


■ 9 4n 








g n3 




gUn 








9 21 


. 


. 








\g 22 


. 


. 





0/ 



(6.6) 



Having established the block structure of the matrix of coefficients of the n— dimensional 
predetermined system g^g^ k = with a homogeneous zero R. H. S., it is now easy to show 
that " ( -" +1 - > equations can be chosen so that the determinant of coefficients will be zero. 



Let us take for example the first 



n(n+l) 



equations from the system of {n 2 — n) equations 



with a zero R. H. S. with the corresponding 



n(n+l) n(n+l) 



dimensional block matrix. In 



the particular case the dimension of the block matrix is determined by the number of block 
matrices on the horizontal ( ) anc [ on ^he vertical ( n ^ n 2 +1 ^ ). Therefore, outside this 
matrix will remain a matrix of " ( " 2 +1 ' ) matrix block columns and (n — 1) — n( " 2 +1 ^ = ( - n ~ 3 ' > 
block rows. The last assumption presumes that the spacetime dimension number n is an 
odd one, so that (n + 1) and (n — 3) are dividable by two. Otherwise, if n is an even 
number, one may consider | x | dimensional elementary block matrices Bk t k ■ Then the 
full block matrix of the system will have (n + 1) block matrices on the block horizontal 
(i.e. (n + 1) block columns) and 2n 2 matrices on the block vertical (i. e. 2n 2 block 
rows). The block matrix of the homogeneous system of equations (with a zero R. H. 
S. ) will be 2(n — 1) x (n + 1) "block" dimensional. The chosen block matrix will be 
(n + 1) x {n + 1) "block" dimensional. Outside this matrix there will remain a block 
matrix of2n — 2 — n — 1 = n — 3 block rows and (n + 1) block columns. 

Let us now compute the determinant of the triangular matrix (5.5), from which we 
take the first ^ (or f ) block rows. This ^ x ^±1) ( or f x f ) block determinant S 
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will be equal to 

n + l 

2 

S =Y[(detB ti ) . (6.7) 

i=i 

But for i ^ 1, 2 the expression for detBu has to be found from formulae (6.3) for the block 
matrix B^. Since only one of the submatrices L\ is different from zero, it is clear that 
detBu = detLi.O = 0, therefore the whole expression (6.7) equals zero. 

Thus we have proved that by taking i^tii (or |) consequent block rows from the 
initial (quadratic and triangular) block matrix, the obtained block matrix will have a zero 
determinant. 



7 DISCUSSION 

The transformation of a matrix system of algebraic equations (with the unknown matrix 
A assumed to be a symmetric one) into a system of linear equations is related to some new 
problems of algebraic nature. The fact that the initial matrix system is a predetermined 
one and solutions are not unique should be reflected in the finally obtained matrix B T B, 
but presently it is not known how the properties of the initial matrix system are "encoded" 
in this matrix. 

It is interesting to mention about some possible physical applications of the proposed 
mathematical algorithm. In theories with extra dimensions, the graviton's field compo- 
nents Iimn(x) in the Kaluza-Klein's expansion under toroidal compactification along the 
extra- dimensional coordinate y are contained the following fields [3]: 

the radion H k = \hf j , (7.1) 
k J 

the scalars s| = h* - ^yyfoy + ^)H k , (7.2) 

the vectors Vj- = > ( 7 - 3 ) 

the gravitons G% = h\ v + |( v + ^)# * . (7.4) 
These fields constitute the (4 + n) x (4 + n) dimensional matrix of the bulk graviton 

(7.5) 



Qk yk 

V k S k 

H3 l 3 



where k is a n— component vector, denoting the Kaluza-Klein (KK) numbers along the 
various extra dimensions. Clearly the matrix (7.5) is a block-structured one and the 
elements in it satisfy the following constraints 

(n + P)G% = ; 0"G£ = O ; of = (7.6) 
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Pvjj = , d»V* = , 



(7.7) 



tiS? k = ,S^ = 0. (7.8) 

The operators (□ + k 2 ) and 5 M may be assumed to have some representation. Then 
the system (7.6-7.8) is similar to the homogeneous part (when % ^ k) of the investigated 
system g^g-* 1 * = 5^. Since the constraints are acting on a part of the elements of the 
block-matrix (7.5), it will be predetermined also. 

Making use of the defining equalities (7.1-7.4) for the radion, the scalars, the vectors 
and the gravitons, the system of equations (7.6-7.8) can be written with respect to the 
graviton's field components Iimn^)- From a physical point of view, a solution in terms 
of these components may be much more valuable. 

The method, however, will require a modification of the approach in the preceeding 
sections, since there it was assumed that the matrix of the coefficient functions A is a 
symmetrical one (g^ = g^), while finding out the graviton's components will be based on 
an initial matrix system, which will be a non-symmetric one. 
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